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Abstract: We investigate up to which order quantum effects can be neglected in cal- 
culating cosmological correlation functions after horizon exit. As a toy model, we study 
(p^ theory on a de Sitter background for a massless minimally coupled scalar field (j). We 
find that for tree level and one loop contributions in the quantum theory, a good classical 
approximation can be constructed, but for higher loop corrections this is in general not ex- 
pected to be possible. The reason is that loop corrections get non-negligible contributions 
from loop momenta with magnitude up to the Hubble scale H, at which scale classical 
physics is not expected to be a good approximation to the quantum theory. An explicit 
calculation of the one loop correction to the two point function, supports the argument 
that contributions from loop momenta of scale H are not negligible. Generalization of the 
arguments for the toy model to derivative interactions and the curvature perturbation leads 
to the conclusion that the leading orders of non-Gaussian effects generated after horizon 
exit, can be approximated quite well by classical methods. Furthermore we compare with 
a theorem by Weinberg. We find that growing loop corrections after horizon exit are not 
excluded, even in single field inflation. 
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1. Introduction 

The precision of measurements of temperature fluctuations in the Cosmic Microwave Back- 
ground radiation has increased enormously in the recent past and is expected to increase 
even more in the near future. From these measurements, statistical properties of the pri- 
mordial cosmological perturbations can be deduced. They are found to have a nearly flat 
power spectrum, and to be close to Gaussian. Non-Gaussian effects (see Q for a review) 
might be detected in the future and can provide a powerful tool to discriminate between 
different inflation models. 

There is therefore a large interest in calculating the statistical properties of the pri- 
mordial cosmological perturbations for different inflation models. In these calculations the 
cosmological perturbations are often parameterized by the curvature perturbation which 
is the perturbation in scalar curvature on time slices of uniform density. This parameteri- 
zation is convenient because of the property that, under certain conditions, perturbations 
in remain constant after their wavelengths have grown larger than the horizon length, i.e. 
after horizon exit. This has been shown for linear perturbation theory in [^, to all orders 
in ( for single field inflation in Q, and nonperturbatively for adiabatic perturbations in 
Q. The latter two references use a derivative expansion and therefore assume that effects 
from wavelengths of the order of the horizon and shorter, are negligible. 

In more complicated models for inflation, e.g. those involving multiple fields, the 
curvature perturbation is not constant after horizon exit. Therefore evolution after 
horizon exit might lead to non-Gaussian effects, which has been investigated in |5|, |6|, |^, |8|, 
|, 0, |n|, 0, [1^, |18|, |0|, ^, H, H, m, H, These investigations 



have been done by solving classical equations of motion, which is assumed to be a good 
approximation to the quantum theory, because quantum effects are presumably negligible 



for wavelengths much longer than the horizon length (see [28| for a recent argument). 

The goal of this paper is to investigate up to which order corrections to cosmological 
correlation functions that are generated after horizon exit, can be calculated reliably using 
classical physics. For this we study (p^ theory on an exact de Sitter background for a 
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massless minimally coupled scalar field (j), as toy model for the curvature perturbation C, 
on an inflationary background. We consider correlation functions with (external) momenta 
much smaller than the Hubble scale H in the quantum theory, using the Closed Time 



Path (CTP) formalism (also known as in- in formalism, see e.g. [29, pQ|), as is also done 



by Weinberg in |31, |3^. Furthermore we formulate a classical theory with statistical 
fluctuations, in such a way that correlation functions in this theory can easily be compared 
with those in the quantum theory. 

In short our results for the (j)^ toy model are that the tree level contributions in the 
quantum theory can be approximated quite well by classical physics (which sounds trivial, 
but we argue that this is not completely so) . However we find that the classical approxima- 
tion is not expected to be good in general for loop corrections in the quantum theory. The 
reason is that loop integrals get contributions from momenta with magnitude up to the 
Hubble scale, and the contributions from loop momenta at scales around H are in general 
not negligible. This is supported by an explicit calculation of the one loop correction to 
the two point function, where we use a small mass as infrared regulator. We argue that 
the classical approximation is not expected to be good at scales around and can there- 
fore in general not reproduce these contributions. An exception is the one loop correction, 
for which we find that the classical approximation can be saved by choosing a suitable 
ultraviolet cutoff. 

We generalize the arguments for the toy model to derivative interactions, and apply 
them to the curvature perturbation C,. This leads to the conclusion that non-Gaussian 
effects generated after horizon exit in multifield inflation models can, up to one loop level, 
be approximated by classical physics. We also compare with a theorem derived by Weinberg 
1 31]. We find that it is not excluded that there are corrections to correlation functions of 
C that grow after horizon exit, even in single field inflation. 

We remark that the problem we are addressing, is related to, but different from the 
problem of the quantum-to-classical transition [33, 34, EBI, p^, 37, 38], that deals with 



the way how quantum fluctuations acquire classical properties by decoherence, and with 
the production of entropy. In this paper we are not investigating how a quantum system 
evolves to a classical system; we are considering a quantum system and a classical system 
separately from each other and investigate how well the classical system can reproduce 
correlation functions of the quantum system. 

In the next section we recall the CTP formalism, using a variation of the Keldysh- 
basis, and apply it to (j)^ theory on a de Sitter background (for other applications of the 
CTP formalism to interacting fields in cosmology see |3^, ||, ||, |5|, ||, 



48, 1^, 51, 52, 53, 54, ^). Subsequently we analyze contributions to correlation 
functions with small external momenta, that are generated after horizon exit (but still 
during inflation). In section ^ we show how correlation functions in classical (ji^ theory 
on a de Sitter background can be calculated perturbatively, starting from given initial 
conditions. The perturbative contributions are graphically represented in a way that is 
similar to the Feynman diagrams of the quantum theory. Next in section Q we argue that 
the classical theory can approximate the contributions from small internal momenta in 
the quantum theory quite well, but that this is in general not the case for large internal 
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momenta. The one loop correction is an exception: we show that by choosing a suitable 
ultraviolet cutoff, the classical approximation can be good. We generalize our arguments 
and conclude in section ^. 

In Appendix |^ details of the quantization are given and a comparison is made with 
finite temperature field theory. Appendix p| contains an argument on closed retarded loops. 
In Appendix |C| diagrams in the classical and quantum theory are compared. To illustrate 
the arguments in this paper, we give in Appendix ^ the detailed calculations of the one loop 
correction to the two point function in the (j)^ toy model, that is generated after horizon 
exit. 



2. Quantum theory 

In this section we set up the quantum theory using the CTP formalism, and analyze the 
contributions to correlation functions that are generated after horizon exit. 
The Lagrangian density of (jy^ theory is 



-g 



(2.1) 



where we are using a metric g^^ with signature — +++. Except when we need the mass 
m as infrared regulator, we take m = and ^ = to obtain a massless minimally coupled 
scalar field. The term 5C contains the counterterms: 



5C 



-5i(j) - ^^zdf^ 



2 ^fn^ 



3!^ 



(2.2) 



We added a linear counterterm 6i to keep (0) = for all times, at one loop level^, hence 
up to order 0{X'^). The potential can be stabilized by adding a cj)'^ term if desired. We use 
a spatial momentum cutoff A as ultraviolet regulator. Then 



4(27r)2 



ln- + 0(A^ 



(2.3) 



where is a renormalization scale. 



2.1 Closed Time Path formalism on a de Sitter background 
2.1.1 CTP formalism 

In a system with time-dependent Hamiltonian H{t), that starts in a state |in) at initial 
time ti, the expectation value of an operator Q at time t > tj is given by 



(Qit)) 



m 



Texp(i / dt'H{t'] 

'tin 



Q 



T exp 



t 

-I I dt' H{t') 

'tin 



m 



(2.4) 



^It can be checked that this is possible by calculating the tadpole diagram, using the F two point function 



of equation (2.24) and the infrared and ultraviolet regulators as discussed later in this paper. 
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where T means a time-ordered product and T an anti-time-ordered product. In the Closed 
Time Path (CTP) formalism (or in-in formalism) this expectation value can also be 
calculated using path integrals, from the generating functional 



Vcp+Vcj)' exp 



; / dt' f d^x (£[(/>+] - £[0-] + + J^cf)-) 



The path integral on the second line can be written in short-hand notation as 



Vcp exp 



dt' j d^x (£[(/>] + J(/)) 



. (2.5) 



(2.6) 



where C is the so-called Schwinger-Keldysh contour which runs from tin to t and back. The 
field <j) and source J are split up in (p'^ , J+ on the first part of this contour, and , J_ 
on the second part, with the condition (t) = 4>~{t)- The integration along the contour 
C explains the name Closed Time Path formalism. The path integral on the first line of 
equation (|2.5| ) imposes that at the initial time tin, the state of the system is given by the 
density matrix p(tin)- Expectation values are then obtained by variation of the sources J+ 
and J_: 



(T {(pixi) . . . (piXn)) T {(j){Xn+l) . . . (piXn+m))) = 

(5"+™Z[J+,J_,p(tin)] 



, (2.7) 



J+,J-=0 



dJ-{xi) ■ ■ ■ SJ-{Xn) 5J+{Xn+l) ■ ■ ■ 5J+{Xn+m) 



where the times x^ are smaller than or equal to the time t used in definition ( |2.5[ ). 

When calculating these correlation functions perturbatively, we need to know the free 
two point functions with all four possible time orderings: 



G-+(x,y)=i(0(x)0(y))W, 
G+-{x,y)=z{^iy)^ix))^^\ 

G++(x, y) = i(T0(x)0(y))(°) = 0(xo - yo)G-+{x, y) + e{yo - xo)G+- (x, y), 
G— {x, y) = 2(T0(x)0(y)) W = 0(xo - yo)G+- (x, y) + 0(yo - xo)G-+(x, y), 

where the superscript (0) denotes the free field correlation functions. They obey the identity 

G++(x,y) + G— (x,y) =G-+(x,y) + G+-(x,y), (2.12) 

and they can be put together in a matrix: 



(2.8) 
(2.9) 
(2.10) 
(2.11) 



G(x,y) 



G++(x,y) G+-(x,y) 
G-+(x,y) G--(x,y) 



(2.13) 



Note that the two point functions depend on the initial conditions via the dependence on 
p{ti) of the generating functional ( [2? 
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In the context of the classical approximation it is useful to transform the (p'^ and (p 
fields to a different basis, which is a variation of the Keldysh basis (see also [^]): 

The Lagrangian density >C — /^[<?5'^] transforms to 

(2.15) 

The free two point functions in this basis can easily be obtained by the transformation 

G.(.. = RG(.. y)K- = ( ^'^^ f ) , (2.16) 

with 

F{x,y) = -'- (G-+(x,y) + G+-(x,y)) , (2.17) 

G^(x,y) = G++(x,y) - G+-(x,y) = ^(xo - yo) (^-+(2;,^) - G+-(x,y)) , (2.18) 
G^{x, y) = G++(x, y) - G^+(x, y) = d{y^ - xq) (G+-(x, y) - G-+(x, y)) , (2.19) 



where we have used identity ( p. 12 ). They obey the equations 



with 



(a + m2 + ^fl(x))F(x,y) =0, (2.20) 
+ m' + mix)) G^'^ix, y) = ^^^i=^, (2.21) 

= -j^= {^/^)g^''{x)d,) . (2.22) 

The and G^ two point functions are often called the retarded and advanced propagators. 
Note that G^(x,y) = G^(y,x). 

2.1.2 Feynman rules on a de Sitter background 

The metric of the de Sitter background is 

ds^ = -dt^ + a^{t) dx^, (2.23) 

where a{t) is the FRW scale factor. The Hubble rate \s H = d/a. In de Sitter space the 
scale factor is a{t) = aQexp{Ht). We will use conformal time r = — dt' /a{t'), which 
runs from —00 to 0. The scale factor in conformal time is a(T) = —1/Ht. 

As initial state p{T{-a) we take the adiabatic or Bunch-Davies vacuum, for Tin — > —00. 
One expects that other choices will give the same results because this state is an attractor 
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state The free field operator with this initial state is given in equation ( [A.5| ) in 

Appendix 1^ and it can be used to derive the free two point functions ( p.l7[ )-( p.l9| ): 

F{k,Ti,T2) = —3 [(1 + k'^TiT2) cosfc(ri - T2) + k{Ti - T2)smk{Ti - T2)] , (2.24) 



2P 

G^{k,Ti,T2) = 6{ti - T2)-p- [(1 + k'^TiT2) siu A;(ti - T2) - k{Ti - T2) cos A;(ti - T2)] 



(2.25) 



and G^{k,Ti,T2) = G^{k,T2,Ti), and where the two point functions depend only on the 
length of the spatial momentum k = |k|. Representing the (/)(^^ field with a full line and the 
0^^) field with a dashed line, the Feynman rules for the two point functions, the vertices 
and the counterterm are^ 

^ ^ =F{k,n,T2), (2.26) 



n T2 


= -iG'^ik,n,T2) 


= -iG^{k,T2,Ti 










= — iA a^(Ti)(5(ri 


- r2)5{Ti - rs). 








T2,' 


= -^-^a^ in) Sin 


- -T2)5{ti - rs), 








Tl ^ T2 


= -ia^{Ti)6{Ti - 





(2.27) 
(2.28) 
(2.29) 
(2.30) 

When a two point function is attached to a vertex, the corresponding time has to be 
integrated over. A closed loop corresponds with an integral over spatial momentum 
/dV(2vr)3. 

2.1.3 Example: equal time two point function 

As an example to which we will return repeatedly, we consider the equal time two point 
function up to one loop level: 

j d?x e''^-''{4>{T, ^)(t){T, 0)). (2.31) 
The tree level contribution is given by 

F{k,T,r). (2.32) 



There is no contribution with the G^ two point function because that vanishes for equal 
times. At one loop level there are contributions 

A B C D 

^-■Q^ --Q-- --C)-- --^ ' 



^These Feynman rules should not be confused with the graphical representation developed in |60| 
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where diagrams A, C and D have also mirror versions which correspond to interchanging 
the endpoints. There is also the diagram 



, (2.34) 

but it vanishes because of the 0-functions in the two point functions. This is an 
example of the general fact that diagrams with no external two point functions vanish, 
as explained in Appendix ^ Furthermore there are diagrams with dashed lines at the 
endpoints, but these diagrams vanish also because of ^-functions. Diagrams with tadpoles 
are canceled by the linear counterterm b\. 

Using the Feynman rules, the diagrams can be translated to integrals. For example 
diagram A translates to 

A 

— --T = {-i?{-i\f r dn a\Ti) r dT2 a\T2) G^(A:, r, n)F{k, r, rs) x 

J Tin J Tin 



/ (^G''^(|k + Pl,n,T2)F(p,ri,r2). (2.35) 



The symmetry factor is 1 because the propagators in the loop are different (diagrams B 
and C have symmetry factor 1/2). The momentum integral is both infrared and ultraviolet 
divergent. The ultraviolet divergence is the same as the one that occurs in Minkowski 
space and is canceled by the counterterm in diagram D. The infrared divergence has to be 
regularized, e.g. by giving cj) a small mass or by taking space to be finite. 

In Appendix ^ the one loop correction ( 2.33| ) is calculated explicitly, using an initial 
time th with \kTH\ < 1- 

2.2 Late times 

2.2.1 Cosmological correlation functions 

In this paper we consider cosmological correlation functions, by which we mean equal time 
correlation functions 

J cfxi . . . (fxr e-^i'i-^i-—*'^'-^'- (</,(t, xi) . . . (P{t, Xr)H-r, 0)), (2.36) 

where the time r is late, i.e. is well after horizon exit with respect to the spatial momenta 
ki, which can be expressed as \kiT\ <C 1. Similarly, early times are times for which \ kiT\ 3> 1. 

We constrain the analysis further by only considering contributions to these cosmolog- 
ical correlation functions, that are generated after horizon exit. For this we introduce a 
split in time at th, a few (Nh) e- folds after horizon exit, such that \kTH\ = exp^—Nn) <^ 1. 
Correlation functions at th have accumulated contributions from earlier times, of which 
we only keep the free field contributions. Then we use these correlation functions as initial 
conditions for the evolution after th- In practice this means that we use the Feynman 
rules as described in the previous subsection, with the only difference that we take th as 



-8- 



initial time. This procedure is not completely correct because the neglected contributions 
generated before th are of the same order in the coupling constant as the contributions 
generated after th- But it does not change the qualitative behaviour of the contribu- 
tions generated after th, and neglecting the contributions from before th simplifies the 
calculations significantly. We comment further on this point in section |5.1| . 

Different contributions to cosmological correlation functions ( p. 36 ) depend in different 



ways on the time r. Those that are proportional to positive powers of r, are dominated 
by their values at the initial time th- These contributions are negligible if th is taken 
sufficiently long after horizon exit. Contributions that are proportional to a non-positive 



power of r can grow after horizon exit and will therefore dominate. In 32 1 it is shown 
that these contributions do not grow faster than powers of ln(— i^r), so negative powers of 
T do not occur. In this paper we call contributions that are proportional to (including 
powers of ln(— ffr)) late time contributions. In this subsection we analyze the dependence 
on T of the different contributions by counting powers of r. 

Contributions can contain integrals over spatial internal (loop) momenta, p, which can 
be arbitrarily large. We have found that the power counting goes quite differently for small 
internal momenta (smaller than the Hubble scale H, i.e. \pt\ <^ 1), than for large internal 
momenta (of the order of the Hubble scale and larger, i.e. |pr| > 1). Therefore we analyze 
first the case that all internal momenta are small, and consider then arbitrary (amputated) 
IPI diagrams for which all the internal momenta are large and the external momenta are 
small. These IPI diagrams can be treated as effective (non-local) couplings in the analysis 
for small internal momenta, and in this way our analysis covers the whole range of internal 
momenta. After this general analysis, we compare our results with the specific case of the 
one loop correction to the two point function. 

2.2.2 Small internal momenta 

When the internal momenta are small, the expressions for the free two point functions F 
and ( |2.24| ), ( |2.25D can be expanded in kn: 



F{k,n,T,)=^[l + 0{k^Tf)], (2.37) 

G''ik,n,T2) =9in-T2)^[k'{rf-Tl) + Oik'T!)], (2.38) 

where k'^rf indicate all possible combinations k'^rf, k'^TiT2 and A;^r|, and similarly for 
k^rf . Using the lowest order of these expansions, it is easy to count the powers of Tj of 
the contribution of an arbitrary Feynman diagram. The F two point function does not 
contribute any factor of r^, and the two point function gives a factor of rf . Furthermore, 
a vertex contributes a factor a^(rj) oc r^"^, and an integral J" drj oc n, so effectively it 
contributes a factor r~^. 

We can divide the Feynman diagrams into two classes: diagrams that contain only 
vertices with one dashed line (diagrams A, B and D in the example of the two point 
function ( |2.33| )), and diagrams that contain one or more vertices with three dashed lines 
(diagram C in ( p.33[) ). Because each dashed line is attached to a G^ two point function. 
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diagrams of the first class have an equal number of vertices as two point functions. 
Each vertex contributes a factor t^^ and each two point function a factor of rf and 
therefore diagrams from the first class are proportional to rf. Because these diagrams can 
contain integrals like J (irj/r,, they can be proportional to powers of ln(— Tj), which are 
largest for the upper limit of the time integrals, i.e. r. 

Diagrams from the second class have more G^ two point functions than vertices, and 
are therefore proportional to positive powers of Tj. For example a diagram with one vertex 
with three dashed lines has two more G^ two point functions than if this vertex would have 
had one dashed line, and is therefore suppressed by a factor of |A;Tj|^. The contributions 
from diagrams of the second class are largest for the lower limit of the time integrals, i.e. 
for th . Therefore the contributions of diagrams from this class are suppressed with respect 
to contributions of the first class by a factor of \kTH\^ = exp(— GA''/^), and they do not have 
the growing factors of ln(r/rj/). 

The power counting of Ti that is done here, is similar to the power counting of the scale 
factor a in the derivation of the theorem in |31]. The expansion of the two point functions 



(2.37) and ( |2.38| ) can be compared with the asymptotic expansions of the wavefunctions for 
late times in |^T| . The difference is that here we use power counting to differentiate between 
growing and vanishing contributions to correlation functions for a specific (c^^) interaction. 



whereas in [31| it was used to differentiate between interactions leading to different late 



time behaviour. We return to this in section 5.2 



2.2.3 Large internal momenta 

We consider an arbitrary amputated IPI diagram, which has small external momenta, and 
where we integrate the internal momenta starting at a scale M somewhat smaller than H. 
We use a cutoff A as ultraviolet regulator. Because both M and A are physical scales and 
not comoving scales, the limits of the momentum integrals are time dependent: Ma(rj) 
and Aa(ri), where Ti corresponds with the time of one of the vertices. We take for this 
time the earliest time that occurs in the loop, because that corresponds with the smallest 
cutoff.^ In this way, the momentum integrals do also contribute time dependencies, which 
we also have to take into account. 

Suppose that our arbitrary amputated IPI diagram has E external lines and vertices, 
each with V legs. Then there are P internal lines and L loops with 

P=^{NV-E), (2.39) 

L = i {NV -E)-N + l. (2.40) 



Furthermore from equations ( p. 24 ) and ( 2.25 ) we see that each internal two point function 



contributes factors proportional to 

(Plln}^e±ip^r,,^ (2.41) 

^If an ultraviolet divergence is local, it occurs only if the times of the vertices in the loop are equal (i.e. 
they are proportional to 5(r,j — Tj^) for all times rj. ); then it does not matter which time one chooses. 
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where = 0,1,2. Each vertex gives / drj/rj, and each loop gives an integral / d?Pi. 
We ignore powers of the external momenta because the internal momenta are much larger, 
k. Then we can count the powers of p and r of the diagram: 



(-3+n)P+3L-l -3N+nP-l 
Pi j 



Aa(Tj) 
Ma{Tj) 



^{NV-Ey3N-l+3 



(2.42) 



where n is the sum of the n^, and where / is a non- negative integer that represents the 
fact that p integrals can also lead to factors of l/r instead of an extra factor p (see e.g. 
the integrals ( p.22| )-( [D.30| ) in the example calculation of Appendix We ignore the 
contribution from the lower limit Ma{Tj), because the full result cannot depend on the 
split in the integrals.^ 

Apparently the power of Tj is independent of the details of the calculation: contribu- 
tions from large internal momenta to IPI diagrams are always proportional to t~^. 

We can now compare the IPX diagram, seen as an effective coupling, with a tree level 
coupling. They are both proportional to t^^, but they differ in the possible numbers of 
(external) dashed lines. The IPI diagrams can have any number of dashed lines, instead of 
one or three dashed lines for the tree level coupling. However, as shown in Appendix]^, it 
turns out that IPI diagrams with no external dashed lines vanish. Hence the non-vanishing 
IPX diagrams can have one or more external dashed lines. 

The IPX diagrams, being effective couplings for small momenta, can be put in the 
analysis for small internal momenta of section p. 2. 2 . The power counting argument of that 
section shows that only those effective couplings with one external dashed line can lead to 
contributions proportional to t^. Effective couplings with more external dashed lines lead 
to contributions that are suppressed by (at least) a factor |A;r//|^. 

The power of A/H in equation ( 2.42D does depend on the details of the calculation. 
When it is non-negative, it can cause an ultraviolet divergence. Some divergent terms are 
proportional to S{Tj-^ ~Tj2) times r,-. that occur in the loop. These divergences are 

the usual local divergences and are canceled by counterterms. Non-local divergent terms 
can also occur (an example of this is given below), as a consequence of the way in which the 
total correction is split up in contributions from individual diagrams. They must cancel 
between different contributions, in order to make the total result finite. 

It is interesting to consider the errors that occur when the ultraviolet regulator A 



is not taken to infinity but kept finite. FVom equation (2.42) it appears that, after the 
counterterms have been taken into account, the errors will be proportional to positive 
powers of H/A. Hence if A is taken to be smaller than H, large errors occur, but if A 
is taken to be larger than H the errors are suppressed. Clearly, internal momenta of the 
order magnitude of the Hubble scale H still contribute to the correlation functions, even 
though H is much larger than the external momenta k/a{T). This is a feature of quantum 
field theory in de Sitter space that is different from what one would expect from field 
theory in flat space. In the latter case, from the point of view of effective field theories, 
one expects only contributions from internal momenta of the order of magnitude of the 



*This is confirmed for the specific case considered in Appendix |d| 



- 11 - 



external momenta. Contributions from higher scales are said to decouple. In de Sitter 
space, scales decouple only when they are larger than the Hubble scale H. 



2.2.4 Example: Late time contributions to equal time two point function at 
one loop 

The general analysis of this subsection can be checked in the example of the late time 
contributions to the one loop correction to the two point function ( 2.33| ), as calculated in 
Appendix 0. In this calculation the momenta are split up between small and large at a 
comoving scale Mem, which obeys |McniT| <C 1 and Mem > k. 



Small internal momenta. The results for small internal momenta of diagrams A and 
B are given in equations (D.ll) and ( p.58| ), respectively. Diagram C does not give any late 
time contribution. After attaching the external lines, the dominant terms can be found in 
equations (D.49) and ( p.71| ) and are proportional to 



e Th 



with 



m 



or 



A^ln^— , 
Th 



(2.43) 



depending on the values of e and ln(r/T//), where e is the infrared regulator. In this 
calculation we have used an expansion that is valid for |eln(— A;r)| < 1, and therefore this 
calculation is only valid for a limited amount of time. When |eln(— A;r)| approaches 1, the 
term on the right in ( |2.43|) becomes of comparable magnitude to the term on the left. 

Two powers of the logarithm ln(r/TH) come from the two time integrals corresponding 
with the two vertices. The extra factor of Iu^t/th) in the term on the right in ( p. 43 ), is the 
consequence of the momentum integration, and was observed earlier in a similar calculation 
in 1 51], where also a small mass was used as infrared regulator. 



Large internal momenta. For large momenta the result for the amputated diagrams are 
given in ( D.42 ) for diagram A, ( D.61 ) and ( D.62| ) for diagram B and ( D.65| ) for diagram C. 
Diagram A has a local ultraviolet divergence that is canceled by the counterterm, diagram 
D. Diagrams B and C have divergent terms that are non-local, and that cancel each other. 
The finite remainder is suppressed for late times. The contribution ( D.62| ) from diagram 
B only removes the dependence on the scale Md. Therefore only diagram A leads to late 
time contributions, in agreement with the result above that late time contributions can 
only come from IPI diagrams with one external dashed line. 

The term from diagram A that grows quickest for large internal momenta, after at- 



taching the external lines, can be found in equation (D.49) and is proportional to 



A^ln^— , 

TH 



(2.44) 



which is comparable to the term on the right in ( 2.43| ) for the small internal momenta. 



Complete result for the one loop correction 

contributions to the one loop correction is two times equation 



The complete result of the late time 
!49| ) added to equation 
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( D.71 ), which gives 



^ f 7 392 7 17 2 4, ^ ^ 2^ 4 
\ 7 TT ln2 - 4Cf3) - In — 



In 2 - 4C(3) - In ^ + - ln(-A;r//)+ 



36(27r)2A;3 \9e 27 3 ' 18 3 ^' H 9 

2 17 208, 2u 8 , , A T 

- + 15 - — 7 - -TT^ - - In 2 - 3 In ^ + - ln(- fcrn) In — + 
e 000 II o / T}j 

- + ^-27 - 2 In 2 + 4 ln(-A:Tj^ )^ ^ + | In^ — + 0( — ) + 0(e) 1 . (2.45) 

e 3 J TH 6 TH TH J 

Note that a consequence of the growing behaviour of loop corrections is that the theory 
become nonperturbative if one waits long enough. 

3. Classical theory 

In this section we consider classical </>^ theory for a massless minimally coupled field (f). 
The evolution of classical fields on a de Sitter background is governed by the equation of 
motion (which can be derived from the Lagrangian density ( |2.1| )) 



dl(t){x) + 2Ha{T)dr(t>{x) - V^cj){x) + a\T)^cl)\x) = 0, (3.1) 

where we use x = (r, x) with r conformal time. Initial conditions have to be imposed at 
an initial time Tin- We focus on the calculation of equal time correlation functions 

((/>(r,xi) . . . (?:)(r,x„))ci, (3.2) 

where the subscript "cl" denotes a correlation function in the classical theory. In this 
section we show how to calculate these correlation functions in a way that is similar to 
the interaction picture in quantum field theory: first we calculate the free field correlation 
functions starting from the initial conditions and using the free equations of motion, and 
then we calculate perturbative corrections, expressed in terms of these free field correlation 
functions. In |^, ^ this method was used in the context of thermal field theory. Fur- 
thermore, we show that the contributions to the correlation functions can be represented 
graphically in a way that is similar to Feynman diagrams. 

3.1 Perturbative calculation of correlation functions 

We assume that at the initial time Tin, initial conditions are given for the correlation 
functions 

(</>(rin,xi) . . . (/)(rin,x„))ci, (3.3) 

and first order time derivatives of these correlation functions. In the free field case (A = 0), 
the initial conditions can be evolved in time using the free field equations of motion. Then 
one obtains the free field correlation functions 

((/)o(ti,xi) . . . (?:)o(t„,x„))ci, (3.4) 

where the subscript "0" denotes the free field solutions, and where the times ri, . . . ,r„ do 
not have to be equal. 
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To calculate perturbative corrections to the correlation functions, we first solve the 
classical equation of motion for (j){x) perturbatively. The first order correction is 

Mx) = -^J d*ya\yo)G''ix,y)cPl{y), (3.5) 

where yo denotes conformal time, and where the retarded propagator G^{x,y) is the solu- 
tion of 

^ {31 + 2Ha(.„)a„ -V-)G«i.,y) = (3,6) 
where we used de Sitter metric ( p. 23 ). This equation is the same as ( ^^ , so that the 



retarded propagator in the classical theory is equal to the one in quantum theory, given in 
equation ( 2.25| ) after a spatial Fourier transform. Higher order perturbative corrections to 



the solution of the equation of motion are obtained by 

Mx) = -jJ d^ya\yo)G''{x,y)Y,Hy)^^~3-l{y)■ (3-7) 
By iteration the i-th order solution can be expressed in terms of the zeroth order solution 



(/)o(x). The full perturbative solution of the equation of motion (3J.) is the sum 



(P{x) = Y,<^i{x). (3.8) 

i 

Perturbative corrections to the correlation function (|3.2[) are obtained by replacing the 



Xj)'s in (3.2) by the perturbative solution (|3.8D, and ordering the terms according to the 



total powers of A: 

(</>(r, xi) . . . (/)(r, x„))ci = ^ (0(r, xi) . . . (/)(t, x„))^i, (3.9) 

r 

with 

((/>(r,xi)...(/)(T,x„))j;i = ^ (0ii(r,xi)...(/)i,,(r,x„))ci. (3.10) 

ii+...+in=r 

When the </'ij's are completely expressed in terms of free field solutions (t)o{x), the cor- 
rections to the correlation function (|3.1[l| ) are expressed in terms of free field correlation 
functions, which we have obtained from the free equations of motion and the initial condi- 
tions in (I 



3.2 Graphical representation 

When we choose the initial conditions to be Gaussian, it is possible to represent the contri- 
butions on the right hand side of equation ( 3.10D graphically in a way that is similar to the 



Feynman diagrams of the quantum theory. In the free field case, Gaussian initial condi- 
tions evolve to Gaussian free field correlation functions. Therefore the free field correlation 
functions are completely determined by the two point function, which we call suggestively 

Fc\{xi,X2) = {(t)Q{xi)(l)Q{x2))c\- (3.11) 
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We assign graphical rules analogously to the quantum case: 



and furthermore 



(3.12) 
(3.13) 

''^ ^A«'(yo), (3.14) 
Mx). (3.15) 



The r-th order corrections on the right-hand side of equation ( p. 101 ) can be constructed 
graphically in two steps. First the (f)iS of equation ( |3.7| ) are represented by tree graphs 
where the endpoint x and the i vertices are connected to each other by propagators. 
The remaining free legs of the vertices are occupied by the 0o's of ( 3.151) . For example the 
second order solution (t)2{x) can be represented by 



.< 02(x) = 2jd^y a\yo) {-i)G''{x,y)^My) x 



d^za\zo) {-i)G%,z)^c^liz), (3.16) 



where the factor 2 comes from two equal contributions. 

In the second step the tree graphs representing the (pi- in ( p.lO| ) are glued together at 
the crosses in all possible ways. When two crosses are glued together, a full line is created 
representing the free two point function ( |3.13| ). Consider for example the contribution 
((/)2(xi)(/)o(x2))ci to the second order two point function {(f){xi)(j){x2))^]^^ ■ The tree graphs 
representing 02(^:^1) and (po{x2) can be glued together in two ways: 

(-i)'(-iA)2 / A a\y^) G«(xi,y) x 



d^za\zQ) G'^{y,z)F,i{y,z)Fa{z,X2), (3.17) 
(where an extra factor 2 comes from two ways of contracting the (po's), and 

O (-^)^(-^A)- 

Xl j X2 

' d^za\zo)G'^{y,z)F,i{z,z). (3.18) 



^2 

d'^y a^ivo) G^{xi,y) Fd{y, X2) x 



The former diagram is equal to diagram A in the quantum theory (equation ( |2.33| ) , after a 
spatial Fourier transform) if Fci{xi,X2) equals the F two point function in equation ( |2.24| ). 
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Both diagrams ( 3.17 ) and ( |3.18| ) can be divergent, depending on Fc\. If the divergences 



are local they can be canceled by adding counterterms. In fact diagram ( |3.18| ) contains a 
tadpole diagram, which is automatically local and can be canceled completely by a linear 
counterterm. 

Similarly one can construct the contribution {(l)o{xi)4'2{x2))dj which is equal to the 
mirror version of diagram A in ( 2.33| ). Finally there is (</'i(xi)(/)i(x2))cii which is equal to 



diagram B in the quantum theory. It is not possible to obtain diagram C in the classical 
theory. 

Note that the resulting classical diagrams can have loops. This illustrates that loop 
corrections occur not only in the quantum theory, but also in the classical theory. These 
diagrams do not vanish, because there are statistical fluctuations. 

4. Classical approximation 

In section ^ we have investigated the late time behaviour of the quantum theory, and in 
section ^ we have set up the classical theory. The graphical representation of the classical 
perturbative corrections as described in section |3^ suggests that the classical theory repro- 
duces exactly the diagrams of the quantum theory with only vertices with one dashed line. 
In Appendix ^ a precise argument is given that shows that this is indeed the case. Hence 
if we choose the initial conditions of the classical approximation such that the classical 
free field two point function F^i is equal to the quantum F two point function ( 2.24| ), the 



classical approximation reproduces the contributions of the quantum theory coming from 
these diagrams. 

In the one loop correction to the two point function ( |2.33D , this means that the classical 
approximation is given by diagrams A, B and D (counterterms are still necessary in the 
classical approximation) . 

In this section we investigate how good the classical theory is as a classical approxi- 
mation to the quantum theory for late times. 

4.1 Small internal momenta 



As argued in section p. 2. 2 for small internal momenta, the diagrams with only vertices with 



one dashed line give exactly the contributions in the quantum theory that are proportional 
to T^. The other diagrams, that have vertices with three dashed lines and are not in the 
classical approximation, give contributions that are suppressed by |A:rj7|^ = exp^—GNn), 
because each vertex with three dashed lines leads to two more retarded propagators com- 
pared to a vertex with one dashed line. These contributions do not grow after horizon exit. 
Therefore, for small internal momenta the classical approximation is good up to errors that 
are suppressed by a factor of exp{—6NH) with respect to the late time contributions. 

This is confirmed in the example of the one loop correction to the two point function. 
Here the late time contributions from small internal momenta are completely coming from 
diagrams A and B and these diagrams do indeed occur in the classical approximation. 

The internal momenta in tree diagrams are always small. Therefore the tree level 
contributions in the quantum theory can well be approximated by a classical approximation. 
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if th is chosen sufficiently long after horizon exit. Note that this is not trivial: the quantum 
theory contains tree diagrams with vertices with three dashed lines that do not occur in 
the classical approximation. 

4.2 Large internal momenta 

In section p. 2.3 we have seen that in the quantum theory, loop corrections get late time 



contributions from internal momenta up to the Hubble scale H, from IPI diagrams with 
one external dashed line. This set of diagrams is not the same as the diagrams of the 
classical approximation. Namely, the IPI diagrams with one external dashed line can 
contain vertices with three dashed lines, but these diagrams do not occur in the classical 
approximation. Therefore, the classical approximation misses late time contributions from 
large internal momenta. This could have been expected: the classical approximation is not 
supposed to be good for physics at scales around H. 

For one loop corrections, the classical approximation can be saved: it turns out that 
the classical approximation does not miss any late time contributions, because at one loop 
level there are no IPI diagrams having both one external dashed line and a vertex with 
three dashed lines. However, another problem arises: because the classical approximation 
has fewer diagrams than the quantum theory, not all ultraviolet divergences are canceled. 
In the next subsection, we treat these new ultraviolet divergences by introducing a cutoff. 

In the two point function, this becomes apparent by the fact that the classical ap- 
proximation does not contain diagram C. As a consequence, the ultraviolet divergence of 
diagram B is not canceled. As mentioned above, this ultraviolet divergence is not local, 
and can therefore not be canceled by a counterterm. 

4.3 Classical approximation at one loop 

In order to deal with the ultraviolet divergences in the classical approximation at one loop, 
we introduce a cutoff at a physical scale A (we use a physical scale and not a comoving 
scale, such that we can use the arguments of section 2.2.3| ). In the following we estimate 



the magnitude of the error that is induced by this cutoff, in one loop diagrams. 

First of all there are the ultraviolet divergences that made the cutoff necessary. From 
equation ( p.42|) we see that they are proportional to (A/iJ)"'^ with ni > (the case ni = 
can give a factor ln{A/H)). As argued above, these divergences occur (at least at one 
loop level) only in IPI diagrams with two or more external dashed lines, and are therefore. 



according to the argument in section 2.2.3 , suppressed by a factor exp{—3NH) or a positive 
power of this factor. Thus errors coming from the divergent terms have an extra factor of 
exp{—3NH){A/ H)'^^ with respect to late time contributions from large internal momenta. 

This can be checked for the example in Appendix 0: the term of the amputated 
version of diagram B that causes a linear divergence in the classical approximation is given 
in ( p. 61 ). Together with the external two point functions, the divergent term has 
an extra factor \kTH\^{A/ H) with respect to non-vanishing late time contributions, e.g. 



equation ( D.62|) 



Secondly there are the errors from terms that are proportional to inverse powers of A 
and that would vanish if the cutoff would be sent to infinity. For IPI diagrams with one 
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external dashed line, these errors can be proportional to and give late time contributions. 
These errors therefore have only an extra factor of {H/ K)^'^ with n2 > 0, with respect to 
other late time contributions from large internal momenta. 
The total error thus has scales like 

with respect to other late time contributions, where q are constants of order 0{\). To 
make this factor considerably smaller than one, th should be chosen long enough after 
horizon exit and the cutoff A should be chosen considerably larger than the Hubble scale 
H\ e.g. if ni = 722 = 1 and ci = C2, A should be chosen of the order of exp{3NH/2)H. 

Instead of introducing the cutoff by hand in the momentum integrals, one can also 
remove the ultraviolet divergences by changing the ultraviolet behaviour of the initial 
conditions. For example one can put a cutoff in the initial conditions: 

I d'xe-'''-^{^{TH,^)ct>{rH,0)),, = ^{l + k^T^H)e{A,^-k), (4.2) 

where Acm is now a cutoff in comoving momentum. The physical scale A that corresponds 
with Acm is |Acm/a(''")|, thus the error estimate (|4.1| ) becomes now 



Another possibility is to choose initial conditions as 



1 



n2 

(4.3) 



where there is not a hard cutoff, but the loop diagrams are made finite (except for the 
tadpole diagram, which does not cause problems because it can be canceled by a local 
counterterm) . The time Tc acts like an inverse (soft) cutoff in comoving momenta, so that 
it should be taken small enough to make t/tc S> 1 and large enough to keep |A;r/f|^(r/Tc)" 
small. 

4.4 Example 

In Appendix |D| we have calculated the one loop correction to the two point function in the 
quantum theory. As an example, we have computed numerically the error that comes from 
using a finite cutoff for a specific set of parameters, as a function of the cutoff. This is 
done for diagram A by taking the upper limit in equation ( D.31| ), adding the counterterm 



(diagram D), and subtracting the term that remains finite in the limit A — > oo (i.e. the 
second line of ( D.42| )). Then the external lines are attached and the times ri and T2 are 



integrated numerically. For diagram B the upper limit in equation ( p. 601 ) is used, and 
there are no counterterms. The results are divided by the full correction ( p. 45 ), and are 



plotted, separately for diagrams A and D, and for diagram B, in figure I. Both errors are 
much smaller than the full correction ( p. 45 ). The error from diagrams A and D is clearly 



decreasing for increasing cutoff. The error from diagram B is suppressed, but increasing 
linearly with the cutoff. These results agree with the arguments given in this section. 
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Figure 1: Plot of the numerically calculated error with respect to the complete result ( 2.45 ), for 
diagrams A and D (full) and diagram B (dashed), versus the cutoff K/H . We have used kru — —0.4, 
kr = -0.03, e = 0.1, and 2^ = H. 



5. Discussion and conclusions 
5.1 Early time contributions 

Both in the quantum theory and in the classical theory we have neglected early time 
contributions, i.e. contributions from times before th- In both cases they can be included 
by imposing initial conditions at th, which can be obtained by calculating equal time 
correlation functions at th in the quantum theory with initial time Tin — > — oo. 

In the quantum theory the initial conditions can be represented by non-local n-point 



vertex functions that act only at the initial time th, as explained for example in |3C]. These 
vertex functions can have any number of dashed and full lines. They can give by themselves 
constant late time contributions, and can also occur in diagrams that give growing late 
time contributions. These extra ingredients make the arguments and calculations more 
complicated, but do not change them qualitatively. For example the calculation of the one 
loop correction to the two point function in Appendix ^ will have extra contributions from 
early times, of order but without growing factors ln(T/r//). 

In the classical theory a similar thing can be done. A practical problem is that not 
all the vertex functions can be represented in the classical theory, because they can have 
any number of dashed lines, while in the classical theory only vertices with one dashed 
line can be represented. Instead the initial conditions at th can be imposed by adjusting 
the free field correlation functions. They are then in general non-Gaussian and cannot 



be represented only by the free field two point function (3.11). As a consequence, the 
calculations become more complicated, but not qualitatively different, as is the case in the 
quantum theory. 

5.2 Generalization to derivative interactions 

Throughout this paper we have used (f)^ theory as a toy model. For other interactions the 
analysis of this paper can be adjusted, which is straightforward for 0" interactions, but 
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less so for derivative interactions, which are of particular interest for cosmology. 

Let us first consider contributions from small internal momenta. A spatial derivative 
leads to an inverse power of the scale factor, or equivalently to a factor of r. As can be 
seen from the power counting argument in section p.2.2| , this extra factor of r suppresses 
late time contributions and prevents that any factors of ln(r/Ti/) occur. 

For a temporal derivative the situation is more complicated. If a time derivative 
dt^ = —Hndri acts on an F two point function, which can be expanded as ( 2.37] ), the 



constant term vanishes, so the result is proportional to rf and as a consequence, late time 
contributions are suppressed. But if a time derivative acts on a two point function, 
which can be expanded as ( p. 38 ), the result is still proportional to r?, and there is no 



suppression of late time contributions yet. Only when there is also a time derivative 
acting on the other time argument , the result becomes proportional to and late time 
contributions are suppressed. 

Hence when an interaction has one time derivative, e.g. (/x/*^, there are still late time 
contributions: if at all vertices (with one dashed line) the time derivatives act at the dashed 
line, no extra factors of Tj appear and late time contributions are not suppressed. But when 
an interaction has two time derivatives, e.g. cp^cj), there is no way to avoid suppression of 
the late time contributions: in a diagram with N vertices (all of this type and with one 
dashed line), there are two point functions and 2A^ time derivatives, so that there 

must be at least one G^ two point function with time derivatives on both sides, or an F 
two point function with a time derivative. 

For large internal momenta, we can reconsider the power counting argument in section 
2.2.3| . A derivative (spatial or temporal) leads in ( 2.42| ) to either an additional factor pr, 



or to no additional factors. Hence derivatives only change the value of the integer I, and do 
not change the arguments of section |2.2.3| . Apparently late time contributions from large 
internal momenta are not necessarily suppressed by derivative interactions. 

In |^T| these matters are treated in a slightly different way. There the time integrals 
are performed first, and after that the momentum integrals, which is a different order than 
employed in this paper. For fixed external and internal momenta, a theorem is derived 
that shows that if the interactions obey certain conditions, the time integrals converge for 
r ^ 0. For the wavefunctions an asymptotic expansion is used, valid for late times (or 
equivalently small momenta). 

Because in this theorem r — > is taken, relative to which all fixed momenta are small, 
and because of the use of the asymptotic expansion, this theorem can be compared with 
our findings above for small internal momenta. They are indeed in agreement with the 
conditions of the theorem. 

5.3 Comparison to stochastic approach 

It is interesting to compare the classical theory as described in section |3|, to the stochastic 



approach ||62| , 63]. In this approach the field (p is also considered to be a classical field with 
statistical fluctuations after horizon exit. The difference with the classical theory of section 
m is that the fluctuations are not imposed as initial conditions at a flxed initial time, but 
are put into the system stochastically at wavenumber H at all times. 
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The stochastic approach has in the recent hterature been used for different purposes, 
e.g. for calculating non-Gaussianities perturbatively in multifield inflation ^, |ll|, 
and for investigating nonperturbative behaviour in de Sitter space that occurs at very 
late times [64, 65], when the factors ln(— r) have grown so large that they overcome the 



suppression by small coupling constants. In the latter case it has been argued |6^, ^ that 
the stochastic approach can reproduce the terms with the largest power of ln(— r) at each 
order in the coupling constant (leading log approximation). 

Below we first discuss the case of only a massless minimally coupled scalar field with 
nonderivative interactions, as in section ^. Then we make some remarks on theories with 
derivative interactions and with other fields than massless minimally coupled scalars. 

5.3.1 Massless minimally coupled scalar with nonderivative interactions 

The stochastic approach does not use a mass as infrared regulator, but uses a finite lower 
limit for the momentum integrals: the classical field (p is defined to contain only modes 
with comoving wavenumber k > H. Physically this corresponds to considering only a finite 
patch of de Sitter space, the size of which increases exponentially by the expansion. So to 
compare with the classical theory of section ^, we need to reformulate the latter using this 
infrared regulator. 

Apart from the different infrared regulator, the stochastic approach makes two addi- 
tional approximations with respect to the classical theory of section ^. First the classical 
field is assumed not to contain modes with wavenumber k > Ha; hence the momentum 
integrals have an upper limit Ha. Second, the wavefunction of the free scalar field ( [A.4| ), 
that is used to characterize the stochastic fluctuations, is approximated by its leading term 
for late times: (t>k,i{T) z ^^^ . In the classical theory of section |3| this is equivalent to 
taking only the leading term of the expansion of the F two point function, as is done in 
equation ( 2.37| ). Also the retarded propagator is approximated by its leading term, as 



in equation ( |2.38| ). 

We now compare the stochastic approach with the classical theory of section |^ for the 
one loop correction of the two point function in cp^ theory, in particular the term with the 
largest power of ln(— r). For the classical theory we can use the calculation in Appendix 
PI we only adjust the infrared regulator. This means that in the calculations for small 
internal momenta we put e ^ 0, and use a lower limit H for the momentum integrals. For 
diagram A this changes equation ( D.ll| ) to 



6{27ryH^{TiT2)'^ V H 

The calculation for large internal momenta is unchanged and gives a finite contribution 
after adding the counterterm of diagram D. Attaching the external lines and performing 
the time integrals gives for the term with the largest power of ln(— r) 

^ ^.ln3^1 (5.2) 



36(27r)2/i:3 { 3 th 
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For the small internal momenta of diagram B we use the integral 

f^'^^ dp fP+^ _ ]_/ lf_ M^^-k 2k 

so that the analog of ( D.58| ) becomes 



Mem 



(5.3) 



Again the integral for the large internal momenta remains unchanged (we use a cutoff as 
described in section |4.3| ) , and attaching the external lines and performing the time integrals 
gives a leading logarithmic term of ln^(r/T/i-). Apparently there is no contribution to the 
In^ t/th term from diagram B when we use this infrared regularization. 

We can calculate the same quantity in the stochastic approach by using stochastic 
sources. However we will not do this, but instead repeat the calculation of above, using the 
additional assumptions of the stochastic approach (i.e. taking the leading order approxi- 
mations for the propagators ( p. 37 ), ( 2.3§| ), and the upper limit k < Ha in the momentum 



integral). We expect that this does not make a difference for the result of the largest 
power of ln(— r) and therefore we interpret the result of this calculation as the result of 
the stochastic approach. 

With the approximations of the stochastic approach, the calculation of diagram A 
reduces to the calculation for small internal momenta done above, with only the upper 
limit changed from Mem to Ha{T2). For the amputated version of diagram A the result is 
then 



, ^o 1 2(r/ -T,^)ln , (5.5) 

and after attaching the external lines and performing the time integrals, the result for the 
largest power of ln(— r) is the same as in our formulation of the classical theory, ( |5.2[ ). For 
the amputated version of diagram B the result is 

-A^ / fc2 , Ha(T2) -k 2k \ , , 



4:{2TT)'^k^H^{TiT2)^\ iJ2 Ha{T2) + k Ha{T2) 

which, similarly to the contribution ( |5.4| ) for our formulation of the classical theory, does 
not lead to \p?t/th terms. 

In this calculation the stochastic approach reproduces the same leading logarithmic 
term (but not the same subleading logarithmic terms) as the classical theory in our formu- 
lation, using a lower momentum limit as infrared regulator. We remark that if we would 
have split the momentum integrals at a physical scale McmO instead of the comoving scale 
Mem in the calculation of the classical theory, all the contributions to the leading logarith- 
mic term would have come from small internal momenta. It is reasonable to expect that 
this remains true for higher orders. Moreover the approximations for the propagators in 
the stochastic approach are then also the same as we made in the calculation for the small 
internal momenta. Therefore we expect that the stochastic approach will give the largest 
logarithmic term at each order in the coupling, consistent with the arguments in [54, 65]. 
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Because of the used approximations, the stochastic approach has fewer problems with 
the ultraviolet than our formulation of the classical theory (see sections |4.2| and 4.3). The 



drawback of using these approximations is that even at one loop order, only the leading 
logarithmic term can be obtained, whereas in our formulation of the classical theory, also 
the subleading logarithmic terms can be obtained at one loop order, as explained in section 



4.3. 



5.3.2 Derivative interactions and other fields 



As argued in section derivative interactions typically lead to positive powers of r and 
therefore to suppression of late time contributions from small internal momenta. There can 
also be interactions with fields that are not massless minimally coupled scalar fields. The 
wavefunctions of these fields are proportional to a positive power of r for late times, instead 
of to T^, as is the case for massless minimally coupled scalar fields. Therefore the leading 
terms of the expansions of the F two point functions of these fields are also proportional 
to a positive power of r. If a diagram contains such an F two point function, this diagram 
cannot lead to late time contributions from small internal momenta. 

However, even if there are no late time contributions from small internal momenta, it 
is still possible that there are late time contributions from large internal momenta. The 
stochastic approach makes approximations that are not valid for large internal momenta. 
Therefore the stochastic approach can have problems with reproducing the largest logarith- 
mic terms correctly in these cases. In |6^] the stochastic approach is applied to the theory 
of a massless minimally coupled scalar field interacting with a massless fermion. Here the 
problem that there are no contributions from small internal momenta is circumvented by 
integrating out the fermion, and considering the effective theory of the scalar field. 

5.4 Application to the curvature perturbation 

The motivation for this work comes from cosmological perturbations generated during a 
period of inflation. As mentioned in the introduction, a suitable parameterization for these 
cosmological perturbations is the curvature perturbation which typically has interac- 
tions involving derivatives. It is of interest to know late time contributions of correlation 
functions of C, and whether these can be approximated by a classical approximation. 

For small internal momenta, one can derive the interaction terms for the various degrees 
of freedom in a specific model of inflation, and then use the conditions of the theorem of 
1 31] to decide whether these interactions can lead to late time contributions. For single field 



inflation (possibly together with M free massless scalar fields) it is shown in |31| that the 
interactions do obey the conditions, and therefore do not give late time contributions to 
all orders. This can be compared with Q, where it is argued, using classical physics, that 
C is conserved after horizon exit. This argument is only valid for adiabatic perturbations, 
which applies to single field inflation, and for small internal momenta, because the large 



internal momenta are removed by a smoothing procedure. Indeed, we found in section iA 
that for small internal momenta the quantum theory can be approximated quite well by 
classical physics. 
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For inflation models involving more fields, there are typically interactions that do not 
obey the conditions of the theorem in |31], and therefore can lead to late time contributions. 
Correspondingly, the perturbations are not adiabatic in these models, hence the argument 
of IQ] does not apply and it is no surprise that (" is not conserved after horizon exit. 

For large internal momenta the situation is different. Since derivatives do not neces- 
sarily suppress late time contributions, it is possible that loop corrections lead to late time 
contributions, even for single field inflation. Moreover, a classical approximation would 
only be able to approximate these contributions up to one loop. An explicit calculation 
should decide on whether these contributions occur or not. The sample calculation in |31| 
seems to indicate that there are no late time contributions, even for large internal momenta, 
but we are concerned about the fact that in this calculation, r — > is taken before the 
momentum integrals are performed. The terms that are discarded in this way, might lead 
to late time contributions. 

Note that the background spacetime in inflation is not exactly de Sitter, but typically 
has a slowly decreasing Hubble constant H. This time dependence should be taken into 
account when deriving the wavefunctions for the fluctuating fields, and when doing the 
time integrals. Moreover, the fluctuations can react back on this background and in this 
way change the time dependence of H. This backreaction can be calculated by considering 



one point functions of the fluctuating fields, as is done for example in pn 53, 



5.5 Conclusions 

We have investigated up to which order corrections to cosmological correlation functions, 
generated after horizon exit, can be calculated reliably using classical physics. We have 
done this by making a detailed study of cp^ theory on a de Sitter background, for a massless 
minimally coupled scalar field (p, as a toy model. 

In the quantum theory we studied late time contributions (generated after horizon exit) 
to equal time correlation functions with external momenta much smaller than the Hubble 
scale. We found that in loop corrections, the loop integrals get contributions from internal 
momenta up to the Hubble scale H. This is different from the intuition from effective field 
theories in flat space, that loop integrals are dominated by internal momenta of the same 
order of magnitude as the external momenta. Our calculation of the one loop correction to 
the two point function supports the argument that the contributions from large internal mo- 
menta (around scale H) are not negligible: they are proportional to ln^(T/T//), whereas 
the contributions from the small internal momenta are also proportional to In^ {t / th) , 
or to In^ {t/th)/€, with e = m^/3H^. 

Furthermore we found that a classical approximation can approximate contributions 
from small internal momenta quite well, but that this does not hold for contributions from 
large internal momenta. This is not surprising, because the classical approximation is only 
supposed to work for physics at momentum scales much below the Hubble scale H. As a 
consequence, the classical approximation is good at tree level, but in general not for loop 
corrections. An exception is the one loop correction, for which the classical approximation 
can be good if an ultraviolet cutoff A > H is introduced. 
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We argued that the results presented for the cf)^ toy model can be extended to derivative 
interactions, and be applied to the curvature perturbation Q. For non-Gaussian effects in 
multifield inflation models this means that at leading order, which is tree level, can be 
approximated quite well using classical physics. Also the one loop corrections can be 
approximated if a suitable cutoff is chosen. But for higher order corrections, the classical 
approximation is not expected to be good. 

For small internal momenta, we found that derivatives tend to suppress late time 
corrections, in a way that agrees with the theorem derived by Weinberg ||3^. However, for 
large internal momenta, late time contributions need not to be suppressed. Consequently, 
it is possible that the curvature perturbation Q is not conserved to all orders after horizon 
exit, even for single field inflation. There might be contributions to correlation functions 
of C that grow after horizon exit, coming from loop corrections. These contributions would 
be suppressed by powers of the coupling constant H/M^\ and possibly also by slow roll 
parameters, but they would be amplified by powers of the number of e-folds In a. 
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Note added 

Simultaneously with this work, [^] appeared on the arXiv, in which the stochastic approach 
is extended to scalar quantum electrodynamics. Similar as in where the fermion is 
integrated out, here the vector field is integrated out and the resulting effective theory for 
the scalar is used for the stochastic approach. 

A. Free scalar field on a de Sitter background 
A.l Scalar wavefunction 

The free field equation of motion for a scalar field on a de Sitter background is 

dl^{x) + 2Ha{T)dr(p{x) - VV(a;) + a^(^) (rn^ + ^R)<P{x) = 0, (A.l) 

where we use x = (r, x) with r conformal time. After a spatial Fourier transformation the 
solutions for the mode functions are (see e.g. p9|) 

^kAr) = -^^Hi^H-kT), (A.2) 
where H^\—kT) are the Hankel functions for a = 1,2, and where v is determined by 

-' = \-^2- (A.3) 
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For massless minimally coupled fields, m = 0, ^ = and we take = 3/2 (we choose v to 
be positive). Then the modefunctions reduce to 

</>MW = ^-£^(l + iA:r)e-*'=- (A.4) 

and (j)k,2{T) = (Pkii'^)- ^'^ ~^ 4'k,iiT) niode function is proportional to 

e~*'^'^/a(T) and is called the positive frequency solution^. Using these mode functions, 
the scalar field operator can be decomposed as 

Hr,^) = / (03 (e^'^-^«k0fc,i(^) + e-*^-''4'A*fe,i(r)) , (A.5) 

where the annihilation operators Qk and creation operators satisfy the commutation 
relations 

[ak, 4'] = (2vr)3<53(k - k'), [ok, Ok'] = 0. (A.6) 

The normalizations are chosen such that 

[(/<(r,x),7r(T,x')] =i<^3(x-x'), (A.7) 

where tt{t,x.) = a'^{T)dr(p{T,x) is the conjugate momentum, and that a(T)(/)(r, x) is a 
conventionally normalized free field for kr — > — oo. The state |0) defined by 

OklO) = (A.8) 

corresponds therefore with the free vacuum state in Minkowski spacetime for kr — > — oo. 
This state is called adiabatic or Bunch-Davies vacuum and is the state with respect to 
which we calculate correlation functions in this paper. 

A. 2 Particle creation 

Despite the confusing nature of the concept of particle number in curved spacetime (see e.g. 
||6^ ), we give here a definition of the particle number in the frame of a comoving observer in 
de Sitter spacetime. This definition of particle number then allows us to make a comparison 
with the classical approximation in thermal field theory in the next subsection. 

Using comoving time, the free Lagrangian is the spatial integral over the Lagrangian 
density (2.1) with A = 

L[ip^,dr^^] = I ^Qa2|5.(^kp-^a2fe2|^kp-^(m2 + ei?)aVk|'), (A.9) 
where we have taken 

/H^k 
(^V'ke^''-, (A.IO) 



^ Often the Hankel functions are used with a negative argument: Hl°'\kT). Then the 4>k,2{T) mode 
function is the positive frequency solution. 
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where the (pi^ are time dependent operators. In this form the Lagrangian describes a system 
of uncoupled harmonic oscillators with frequencies ujk = \/k'^ + a?{m? + ^R). The conju- 
gate momentum is defined as TTk = 5L/6{dr'P]i) = a^dr^k, and using this the Hamiltonian 
becomes 

d'^k /^l|vrkP ^ 1 2t2| |2 I ^/ 2 I cD\ 4| |2 



H[^k,nk]= J ^i^-L-^ + -a'k'\ipk\' + -{m'+^R)a^\ipk\'j. (A.ll) 

By analogy to the harmonic oscillator, we define time dependent annihilation and creation 
operators ak and by 



a^k = ^ (ak e-^^ + al^ e"=^) , ^ = i^^f Sk e^^^ - e^'^^ ) . (A.12) 



They obey the commutation relation 

[Qk,4'] = (2^)''5'(k-k'). (A.13) 
The Hamiltonian ( A.ll ) expressed in terms of these operators has the familiar form 

H = J ^;^(^alak + ^[ak,al]^uJk. (A.14) 



If we take m = and ^ = we can use the modefunction ^fc,i(T) of equation (A. 4) to 



express the operators Ok, aj^ ™ terms of the operators aki Ok mode decomposition 



— i g2ifcr 

"k = + 2iA;r)ak + i^^aL,^, (A.15) 

«U = + ^(1 - 2^/cr)al^. (A.16) 

For /cr — > — oo the cik operator becomes equal to ak- We can define and by 

(0|4ak'|0) = nfe (27r)35(k - k'), = (A.17) 

(0|aka-k'|0) = (nfe + ^)(2^)35(k-k'), '''^ = 1^2 ^^-^^^ 

where one can interpret as the particle number and kind of off-diagonal particle 

number, with respect to the vacuum at kr — > — oo. Clearly these quantities are proportional 
to a^. 

A. 3 Comparison with thermal field theory 

It is interesting to compare with thermal field theory on a Minkowski background (see 



1 61, p%). The F two point function in a thermal system with temperature T in Minkowski 



spacetime is given by 



F(/c,ti,t2) = ^(^nth(A;) + 0cosA;(ti-t2) ^th(fc) = ^^/t _ (^-1^) 
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where nth(fc) is the particle number. For small momenta with respect to the temperature 
k <^ T, TT-th — T/k, which becomes large and classical behaviour emerges. Moreover the 
two point function does not have this amplification factor; it is given by 

G'^{k,h,t2) = smk{h-t2). (A.20) 

k 

Therefore a diagram containing a vertex with more than one dashed line is suppressed with 
respect to the same diagram with a vertex with one dashed line. 

To compare this with the de Sitter case, it is natural to consider a0fc(r). Then we have 



a{Ti)a{T2)F{k,Ti,T2) = ^ 



1 1\ \ Tl ~ T2 . . , < 

+ - COS k[Ti - T2) + — sm k[Tl - T2] 



2k'^TlT2 2) ' ' 2kTlT2 

(A.21) 

We see that for equal times ti = T2 = t this quantity grows as 

111 ^ 
+ -ocn + n + -. (A.22) 



2k^T^ 2 2' 

The quantity a{Ti)a{T2)G^{k,Ti,T2) does not have this amplification factor for late times, 
as can be seen from expansion for small kr. Therefore a diagram containing a vertex with 
more than one dashed line is suppressed with respect to the same diagram with a vertex 
with one dashed line, as in thermal field theory. 

Both in the de Sitter case as in thermal field theory, the arguments given above explain 
why the classical approximation is good for small physical internal momenta, i.e. k/a <^ H 
(|A;t| ^ 1) or /c ^ T. But, as we argue in this paper for the de Sitter case, for large internal 
momenta (~ or « T) problems arise for the classical approximation, e.g. in the thermal 



case involving Hard Thermal Loops |70| 



B. Amputated diagrams with no external dashed hnes 

The contribution of an amputated diagram with no external dashed lines, as for example 



occurs in diagram ( 2.34 ), vanishes. The reason is that in such a diagram there is always 
a closed loop of G^ two point functions. This can be seen by picking an arbitrary vertex, 
and from there following one of the dashed lines, along the G^ two point functions, from 
this vertex to a next vertex. From this new vertex one can repeat this to go to the next 
vertex. Because there is no external dashed line, this can be repeated indefinitely while 
staying inside the diagram. Since there is only a finite number of vertices in the diagram, 
one eventually ends up at a vertex for the second time. Therefore there must be a closed 
loop of G^ two point functions in a diagram with no external dashed lines. 

Because of the ^-functions of the G^ two point functions, the times of the vertices of a 
closed loop of G^ two point functions have to be ordered. But in a closed loop of there is 
always at least one G^ two point function for which the ^-function vanishes, and therefore 
the complete diagram vanishes. Even if all the internal times of the loop are equal, the 
diagram vanishes because the G^ two point function vanishes if the time arguments are 
equal. 
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C. Correspondence between diagrams in quantum and classical theory 



In this Appendix we show that, if one chooses initial conditions such that 

Fei(k,ri,r2) = F(k,ri,r2), (C.l) 

the classical approximation reproduces the diagrams of the quantum theory with only 
vertices with one dashed line. 

Discarding the vertex with three dashed lines in the quantum theory corresponds to 
discarding the term with [(p^'^^)^ in equation (2.15). An arbitrary equal time correlation 



function, calculated up to order n in the coupling A, can then be written as 

(<A(r,xi) . . . 0(r,x,)) = (.^(i)(t,xi) . . . 'A^'^(t,x,)1 (Sint)")^"), (C.2) 

where the superscript (0) denotes the free correlation function, and where the modified 
interaction part of the action is given by 

Ant = -f dra^r) J d\ A (</,(!) )2<^(2). (C.3) 

On the right hand side of equation ( p. 21) the fields c/)^^^ and (/»^^-' have to be contracted in 
all possible ways: every (/>(^) is contracted with a </>(^) to a retarded propagator G^, and 
the remaining i;A^^)'s are contracted with each other to F two point functions. Suppose 
that we do have contracted the (/)^^^'s, but not yet the remaining (/)(^^'s. The correlation 
function can then be drawn as a number of disconnected diagrams, in which the vertices are 
connected by two point functions and which have loose ends representing (p'^^^^s that are 
yet to be connected. In total there are n vertices. Each disconnected diagram can be a tree 
diagram or a loop diagram. A loop diagram with only G^ two point functions vanishes, 
as explained in Appendix |^. Hence if one of the disconnected diagrams contains a loop, 
the complete contraction does not contribute to the correlation function, and for all non- 
vanishing contractions the disconnected diagrams are tree diagrams. These tree diagrams 
are the same tree diagrams in the classical theory that represent the perturbative solutions 
(pi in the classical theory ( |3.7| ), where all the (p^^^^s in the quantum theory correspond with 
free field solutions (pQ in the classical theory. The symmetry factors are also equal because 
they arise in the same way (the factor 1/n! is canceled by a factor n\ from the n equivalent 
vertices). The remaining contractions of the (/>^^^'s in the quantum theory are equal to the 
contractions of the 0o's in the classical theory, because of equation ( p.l| ). 

Therefore the classical theory with the same couplings as the quantum theory, and 



using initial conditions such that (C.l) holds, gives all the diagrams using only the vertex 
with one dashed line in the quantum theory, up to vanishing diagrams. Hence this classical 
theory reproduces the late time contributions for small internal momenta. 

D. One loop correction to two point function 

In this appendix we calculate the one loop correction to the two point function as given in 
equation ( p. 33 ). We first consider the diagrams with one external dashed line (diagrams 
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A and D), and then the ones with two external dashed Unes (diagrams B and C). The 
complete result is given in equation ( p.45 ). 



D.l Diagrams A and D 

We start with the diagrams with one external two point function. First we calculate 
the amputated diagrams, and then attach the external lines. The amputated diagrams are: 



k ^ A k D 

) -"^■® (D.l) 

p 

The amputated version of diagram A is given by 

Aamp(fc,ri,r2)= ~^^3~;y j ^^^6H^-p-p')G''{p',n,T2)F{p,n,T2) = 

/ dpp / dp'p'G'^{p',Ti,T2)Fip,n,T2), (D.2) 

Jo J\v-k\ 



where we have used the identity 

d^pd'p'5'{k + p + p')f{k,p,p') = — dpp dp'p'f{k,p,p'). (D.3) 

\p 

We will evaluate this integral below. For diagram D we see from equations ( |2.3[1| ) and (^ 
that it is equal to 

—iX^ /A 

Damp(A;,ri,r2) = -ia^(ri)5„5(ri - T2) = ^^^^^2^4,^4 [j; ) '^('^i " '^2), (D-4) 



where A is the ultraviolet momentum cutoff and /i is a renormalization scale. The coun- 
terterm 6z is finite and leads to terms proportional to positive powers of r, and is therefore 
left out. 



We calculate the integral (|D.2| ) by splitting the p integral in a small momentum part 
jjVfcm ^ large momentum part dp, with |McniTi| <C 1 and Mem > k.^ 

D.1.1 Amputated diagram for small internal momenta 

The integral in equation ( p.2| ) is infrared divergent for p — > 0. We regulate this divergence 
by giving the field a small mass m <^ H, such that u = 3/2 — e with e = m^/3i7^. The F 
and two point functions are then, using equations ( p. IT] ), ( |2.18| ) and 



F{k,n,r2) = ^ . Re(gW(-fcn)gW*(-fcr2)), (D.5) 

4a{Ti)a{T2) V / 

G^(A;,n,r2) = -^^|^^0(Ti-T2)Im(i7«(-A;ri)i/«*(-A;r2)). (D.6) 



''In a similar split of integrals is used to calculate a similar integral. Note however that the integral 
there differs from the integral here, because the self-energy kernel of is not the same as the amputated 
diagram A. 
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Using (see [71|) 

Hl^\-kT) = M-kr) + i (^^^J^{-kT) - -^J^,{-kT)) , (D.7) 

\ Sm VTT Sm UTT J 

and the identity r(z/)r(l — u) = t: / sinz^vr, we obtain 

F{k,Tl,T2) = ^{k\lT2)\ (D.9) 

G^{k,n,r2) = ein - r,)^ [rl (^)^ - . (D.IO) 

The integral is 

' ' dp— j dp p = 



{2T,YkmTt4 6 Jo J\p-k\ 



2t 

cm 



+ 2rf In iMcmTsI - 2r| In |Mcmri| + 0(e) ) . (D.ll) 



Q{2'KfH^{TiT2Y 

D.1.2 Amputated diagrams for large internal momenta 

For large momenta we approximate the field to be massless and we use the two point 
functions of equations ( p.24| ) and ( p.25 ), which we write as 

F(/c,ri,r2) = — ^F,(A;,ri,r2), (D.12) 

i=l 

Fi{k,Ti,T2) = ^ COS kAr, 

F2{k,Ti,T2) = -^ATsinkAr, 
1 

F3{k,Ti,T2) = -TiT2COskAT, 

with At = ti — T2, and similarly 

3 

G^(A:,ri,r2) =0(Ti-r2)i?2^Gf(A:,ri,r2), (D.13) 

i=l 

Gi{k,Ti,T2) = psin/cAr, 
G^(/c,ri,r2) = -^At cos kAr, 

G§{k,Tl,T2) = ^TiT2SmkAT. 
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In the following calculations we use the definitions 

si(x) = r dx' ^ , ci{x)=- r dx' ^ , 

Jo ^ Jx ^ 



which behave for small respectively large arguments as 



Si(x) =x + 0{x^), 



Ci{x) = 7 + lnx- — + 0{x^), 



Si(x) = 2 
Ci{x) = 



TT cos x sm x 



sm X cos X 



+ 0{x-^), 



(D.14) 

(D.15) 
(D.16) 



and the identities 

P'^'' , , sinp'Ar 
dp 



rp+k 
Jp—k 



p 



n 



sin(p + A;) At ^ sin(p — /c)Ar^ 



rp+k 
•Jp—k 



p + k p — k 

Ar(^Ci((p + A;) At) - Ci{{p - A;)Ar)) , 



, , cos y At 
dp 



P 



/2 



cos(p + A) At ^ cos(p — fc)AT^ 



p + k p — k 

- At (Si{{p + A;) At) - Si((p - A) At)) , 



p+k _i 2 

dp sinp At =— — (cos(p + A) At — cos(p — A) At) = -— sin A; At sin p At, 

p-k 

p+k J 2 

dp' cos At =-r— (sin(p + A;)At — sm(p — k)AT) = — — sin AAtcos^At. 

p-k 



At 



Next we calculate the contributions 

rAa{T2) 
'M, 

for i and j from 1 to 3: 

#1: Fiip,n,T2)G^{p',n,T2) 



rAa{T2) rp+k 
/ dppFi{p,Ti,T2) / dp p Gf{p,Ti,T2 

J Mem J p-k 



(D.17) 

(D.18) 
(D.19) 

(D.20) 

(D.21) 



Aa(r2) cospAT /•^'+*^ , , sinp' At 
dp 5 — / dp 



Jm,^ P^ Jp-k P'^ 



cos pAr 
P 



sm{p + k)AT ^ sm{p — A)At^ 



p + k p — k 

AT(ci{{p + A:) At) - Ci{{p - A) At)) 



Aa(r2) 



+ 



-I Mc, 



/•Aa(T2) C 
/ dpi 

J Mom I 



^ sinpAT ^ sin(p + A;) At ^ sin(p — A;)At^ 



P 



p + k p — k 

AT(cii{p + k)AT) - Ci((p- A;)At))^ + 

cospAT / sin(p + A) At sin(p — A)At\ 
P V (^^ + (p ~ J 



(D.22) 
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#2: Fi(p,ri,T2)Gf(p',ri,r2) 

^^"(•^2) cog^^^ fP+k COS p' At 

-k 



/•Mt2; cos p At , , 

-At dp ^ — / dp' 



At 



p- 

cos p At 



P 



P 

/•Aa(T2) 

At / dp 

JMc, 



Aa(T2) 



Ci((p + k)AT) - Ci{{p - A;)Ar)) 
r ^smpAr j^^.^^^ ^ ^^^^^ _ ^.^^^ _ ^^^^^^ 



+ 



cospAr / cos(p + A;) At cos{p — k)AT 



P 



p + k 



p — k 



(D.23) 



#3: Fi(p,ri,r2)Gf(p',ri,r2) 

'■^'^(^2) cos p At , 



T1T2 



dp ■ 



/ dp' sin At = 

Jp—k 



T1T2 sin A;At 

sin 2pAT 
dp ■ 



At 



J Ma 



P^ 

T1T2 sin A; At 
A7 



sin 2pAT 
P 



+ 2At Ci(2pAT) 



Aa(T2) 



(D.24) 



#4: F2(p,Tl,T2)Gf(y,Tl,T2) 

/■A«(r2) sinpAT fP+'' , .sinp'AT 
At / dp — - — / dp' 



Mc, 



/ 



p 

Aa(T2) 



p—k 



p 



,/2 



c/p At- 



sinpAT / sin(p + fc)AT sin(p — A;)At 



p + k 



+ 



p—k 



At(c\{{p + A:) At) - Ci((p - A)AT)j 



(D.25) 



#5: F2ip,Ti,T2)G§{p',n,T2 

J 

J A 



Aa{T2) sinpAT fP+^ , .cosp'At 
At^ / dp / dp' 

Men, P Jp-k P 



Ar 



/•At 

J Ma 



Aa(T2) 



dp 



sin p At 
P 



(Ci{ip + k)AT) - Ci{{p - A;)At)) , (D.26) 



#6: F2(p,Tl,T2)Gf(p',Tl,T2) 



/.Aa(r2) sin p At , . ,^ 

Att\T2 / dp / dp sinpAT = 

^Mcm P Jp-k 

_ /-AaCra) 1 _cos2pAT 

T1T2 sm kAt / dp = 

Jm^^ p 



T1T2 sin A; At 



Inp - Ci(2pAT) 



- Aa(T2) 

- Mem 



(D.27) 
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#7: F3{p,Ti,T2)G^{p',n,T2 
/•A( 

: / dp cos p At / dp' 

-'Mem 

l'Aa(T2) 

dp cos p At 



T1T2 



Aa(T2) rp+k sin At 

Jp 

Aa(T2) 



jp-k 



riT2 



Mc 



sin(p + /c)Ar ^ sin(j» — A;)At^ 



T1T2 



p + k p — k 

Ar(ci((p + A;) At) - Ci{{p - /c)At)) 

sin/cAr(^ln(p2 - A;^) - Ci(2(p + fc)Ar) - Ci(2(p - A:)Ar)) + 
cos fcAr (Si(2(p + A;) At) - Si(2(p - fe) Ar)) 
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Tir2Ar / dp cos p At [Ci{{p + k) At) -Ci{{p- k) At)), (D.28) 

-/Mem ^ ^ 



Aa(T2) 



+ 



J Me, 



#8: Fs{p,Ti,T2)G^{p',Ti,T2) 



Ma(r2) .p+fc COS y At 

riT2Ar / dp cospAr / dp 

J Mem Jp-k P 

fAa(T2) 



T1T2AT dpcospAT[Cii{p + k)AT)-Ci{{p-k)AT)), (D.29) 

-^Mem ^ ^ 



#9: F3(p,ri,r2)Gf(p',ri,r2) 



/■l\a[T2 ) pp+K 

r^r| / dp cos p At / dp' sinp'Ar 

J Mem Jp-k 



^2^2 

^ ^ sinfcAr 



Ar2 



sin^ pAr 



Aa(T2) 



. (D.30) 



Me, 



Together this becomes 
cospAr / sin(p + fe)Ar sin(p — fc)Ar 



rir2 sin /cAr 
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p — k 
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For the upper hmit the boundary term vanishes as 1/A^, except the last term which we will 
discuss below. The lower limit of the boundary term gives (where we use that iMcmTjl <^ 1) 

2. 

3' 



^A;At^ - 2A;tiT2 At( - 2 + 7 + ln2McmAT) + 0(t/). (D.32) 
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Using Mathematica, the integral in ( p.31 ) becomes for |McmTj| <^ 1 



8 - 67 - 61n2McmAr) + 0{4). 



(D.33) 



Together equation ( p.31|) becomes 



2k 



7 



^ (rf-r|)(^--7-ln2Mem(ri-r2)j--/crir2(ri-r2)+^sinA;AT sin^ Aa(r2)Ar+0(rf ). 

(D.34) 

The term that contains the sin^ A is logarithmicahy divergent for A ^ 00. This can 
be seen as follows. Consider the integral 



dAT0(AT)/(AT 



sin^Aa(r2)AT 1 



1 — cos 



dAT/(Ar)- 



-2A At 
H ti-At 







At 



f 

where /(Ar) is a test function. The integral can be split up into two integrals 

lim 

£^0 







, (D.35) 
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where r] is used as a regulator time, which we take to zero in the end, after taking the limit 
A — > 00. In the first integral we can approximate 



Ar 



Ti — At Ti 



— , /(Ar)p^/(0), 
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so that it becomes 



1 — cos 



lim / dArfiO)- 



( -2A Ar 



At 



/(O) 7 + ln- 
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where we have taken Ci{-2Ar]/ Hti) 0, and Ci(-2Ae/i?ri) ^ 7 + ln(-2Ae/FTi). The 
remaining integral is 

l-cos(^^ 



lim / dAr/(Ar)- 



H ti~-At' 



At 



Ar 



(D.39) 



where the term with the cosine vanishes, provided that the test function /(Ar) vanishes 
sufficiently fast as Ar — > 00. Together we obtain for A ^ cxd 



dATe{AT) /{At 



sin^ Aa(r2)AT 

' a7 

dAT/(Ar)- 
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which is in the language of distributions 
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+ 6{At) 7 + ln 
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Using this result in equation (|D.34 ), gathering the right prefactors and adding the 
contribution from the counterterm ( p.4D , we obtain for the large momentum contribution 



^^'^y\~ ^'\, ( -{rf - rl) f- - 7 - ln2M,^Ar) - -nr2(n - r2) + 
2{2tt)^H'^{tiT2)'^\3^ ^ ^^\3 ' cm y ^ i a i 2)^ 



(D.42) 



D.1.3 Attaching the external Unes 

Adding the small and large momenta contributions, we obtain for the amputated diagrams 
A and D: 



{k,Ti,T2) + Damp(fc,n,T2 



2riT2(ri -r2) +2Tf In 
'0{-rj + Ti -T2 



T2 



2(ri - T2) 



2t^ In 



2(ri - T2) 



+ 



Tl - T2 



+ 6{ti - T2) (^7 + In -J^) 



+ 0{T^) + 0{e)\, (D.43) 



where the dependence on Mem has dropped out. The full correlation function is obtained 
by 

-i / dri / dT2G^{k,T,Ti)F{k (/C,ri,r2) +Damp(^,Ti,7"2))- (D.44) 



Because the external momentum k is small, i.e. \kTi\ ^ 1, we can use the expanded versions 
of the two point functions ( |2.37| ), ( p.38|) (or the ones of (^ ), ( plO|) , but this gives only 
corrections of order 0(e)). Using the integrals 



(r^ - rf)(ri - T2) 



dri / dv. 
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+ 6{ti - T2) (^7 + In j 
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(recall that ry is sent to zero), this becomes 

r 1 194 7 17 2 2 , , 1 2u 

■<! - + ^ - 7T7 - T^TTvr^ - -ln2 - 2C(3) - -ln-^ + 



36(27r)2A;3 1 3e 27 6 ' 36 3 ^ 2 H 

2 , 127 17_ 1_2 4,_^^ 3,_^ 2/i^^ r 



H 7 vr^ In 2 In — In h 

3e 18 6 ' 3 3 2 H J th 

^ + ^ - 7 - In 2) In^ — + i In^ — + 0(— ) + 0(e) I . (D.49) 

There is an equal contribution from the diagram with ri and T2 interchanged. Note that 
there is no dependence on In k/fi for |A;r| ^ 1. 

D.2 Diagrams B and C 

The amputated versions of the diagrams with two external propagators are 

B C 
)■- ;-- . (D.50) 

n )^ T2 n \^ T2 

They translate to 

Bamp(fc,n,r2) = 2ff8(^^!^^)4 / -^^l" " P " P')^(p'> ^ , r2)F(p, Ti , Tg) 

\2 /"oo rp+k 

/ dpp / dp'p'F{p',n,T2)F{p,Ti,T2), (D.51) 
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S(2'Kfkm{TiT2Y 



(D.52) 



where both diagrams have a factor 1/2 for symmetry. Diagram C has an additional factor 
1/4 from the vertex with three dashed lines ( ^.29 ). We split the p integral again in a small 
momentum part and a large momentum part. 

D.2.1 Amputated diagrams for small internal momenta 

For small internal momenta we use the expanded propagators ( D.9| ) and ( D.1C| ). 

Diagram B. The integral is 

HHrir2r f'''"^ dp T'dp' ^^^"^'^ - 



2{2^f kH\T^T2Y 4 7o J\p-k\ {PP'Y 



{TlT2y 



8{2TT)^kH^{TiT2)^{2e - I] 



( I dpp-'+'^ {{p + -{k- + 

dpp-^+^' {{p + k)-^+^' -{p- ) . (D.53) 



-37- 



The integral on the middle line of ( D.53| ) is finite, but the individual parts are infrared 
divergent. Therefore we calculate the individual parts for e > 1/2, and in the end use 
analytic continuation to e <C 1. The integrals are (using p = kx) 



dxx-^+^' (l + x) 



-l+2e 



2e 



+ ln2 + 0(e), 



1 
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where we have used analytic continuation in the first two integrals and y = 1/x and z = 1—y 
in the latter two. The right hand side of equation ( p. 53 ) becomes 



2e 



-A2 ^l 



-2 

— + 4 
e 



2v^ + ln, -A_£Hl + o(e; 
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1 , Mem + ^ 
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+ 2 1n(fcVir2) + 0(e) . (D.58) 



Diagram C. From equations (D.52) and ( p. 10 ) we see directly that diagram C does not 
give late time contributions and also does not have an infrared divergence. 



D.2.2 Amputated diagrams for large internal momenta 

The contributions from large internal momenta can be calculated in a similar way as is 



used for diagram A in section D.1.2. 



Diagram B. The sum of integrals 
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(where rg = ti,T2, depending on which time is earher) , is equal to 
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The only ultraviolet term comes from the last term of the boundary term and is, including 
the correct prefactor: 



-A^ sink At 



p + 



sin 2pAr 
2Ar 



Aa(-r/3) 

Mem 
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The only term that gives late time contributions is the first line in the integral. It is 
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Only the last term of the boundary term is ultraviolet divergent: 
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(D.65) 



J Mc: 



The diagram with the vertices exchanged gives the same result, except that 9(ti — T2) is 
replaced by ^(r2 — ri). There are no further late time contributions. 
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Ultraviolet divergences. The ultraviolet divergent terms of diagrams 

(D.66) 



(D.65), and C with the vertices exchanged, add up to 

— sink At sin2pAr 



8{27rykH^{TiT2)^AT At 
which is finite and does not give late time contributions 
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D.2.3 Attaching the external lines 

Adding the small and large momenta contributions, we obtain for the amputated diagrams 
B and C: 

Bamp(fc,n,7"2) + Camp(fc,ri,r2) 



-A2 



4(2.) W W V e + ' ^"^''"^"^^ + ""^"^^ + ""^^V ' 

(D.67) 

where the dependence on Mem has dropped out. The full correlation function is obtained 

by 

-/ dri / dT2G^{k,T,n)G^{k,T,T2)(B.,mp{k,n,T2)+C^mp{k, Ti,T2)). (D.68) 

Jtu JTfl 



Because the external momentum k is small, i.e. |/crj| <C 1, we can use the expanded version 
of the propagator ( p.38| ) (or the one of ( D.IO ), but this gives only corrections of order 
0{e)). Using the integrals 



dr2 ^ T^^^T^ ^ = - + - In — + ln2 — + 0( — ), 
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this becomes 
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